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We give a geometric formulation of the Fokker-Planck-Kramer equations for 
a particle moving on a Lie algebra under the influence of a dissipative and a ran- 
dom force. Special cases of interest are fluid mechanics, the Stochastic Loewner 
Equation and the rigid body. We find that the Boltzmann distribution, although 
a static solution, is not normalizable when the algebra is not unimodular. This 
is because the invariant measure of integration in momentum space is not the 
standard one. We solve the special case of the upper half-plane (hyperboloid) 
explicitly: there is another equilibrium solution to the Fokker-Planck equation, 
which is integrable. It breaks rotation invariance; moreover, the most likely 
value for velocity is not zero. 

1 Introduction 

Arnold [Q showed that the Euler equations of an ideal fluid describe geodesies 
in the group of volume preserving co-ordinate transformations of the space filled 
by the fluid. The negative curvature of this geometry implies that the geodesies 
diverge from each other exponentially: an elegant explanation for the observed 
instability of fluid motion. A more realistic description would include dissipation 
(Navier-Stokes) as well as small random forces not included in the ideal model. 
Since the infinite dimensional fluid system is difficult to understand directly, we 
look for analogous finite dimensional examples that have (i) Negative Curvature 
(exponentially diverging trajectories), (ii) Dissipation and (iii) Random Forcing. 

The celebrated Langevin equation for the velocity of a particle in Brownian 
motion includes (2) and (3); being one dimensional, it cannot have curvature. 

dv 

Here 7 is the dissipation constant, proportional to the viscosity of the medium 
in which the particle is moving. Also, r\ is a random force, modelled as a Gaus- 
sian of zero mean and zero correlation time("white noise"): 



<r](t)r)(t') >= 2D5(t-t'). 



The constant D measures the strength of the fluctuations. This leads to the 
Fokker-Planck equation for the probability density of velocity 



dP _ d 
dt dv 



dP 

D — h jvP. 

ov 



The equilibrium (static) solution of this equation is the Boltzmann distribu- 
tion 

P(v) = \e~^ 
with an inverse temperature given by the Einstein relation 

f3D = r 

How would things change if the particle is moving in a space of negative 
curvature? We will study in this paper the simplest such model. We will show, 
by explicit solution, that the equilibrium solution is not always the Boltzmann 
distribution. In particular, the most likely value of velocity is not zero. The 
essential reason is that the invariant measure on the space of velocities is not 
translation invariant: the Boltzmann distribution is not integrable in this mea- 
sure. There is another static solution to the Fokker-Planck which is integrable. 
In simple cases we can find it explicitly in terms of classical functions. 

There is already an extensive literature on Brownian motion on curved man- 
ifolds. The main approach is to solve the heat equation in position space. How- 
ever, this does not address the question of the equilibrium velocity distribution, 
which we study here. 



2 The Fokker-Planck-Kramers Equation 

We review theory of a dynamical system subject to dissipation and fluctuation, 
paying special attention to the role of the measure on phase space in the Fokker- 
Plank-Kramers equation. We base our discussion on refs. [2j |3] but use a more 
geometric formulation. The erudite reader might want to skip this section after 
a glance at the notation. 



2.1 The Ideal System 

Consider a dynamical system with equations of motion 

du a _ 
dt ~ Va 

for some vector field V a . We will assume that there is a measure [i that is 
invariant under this evolution, 
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d a Wa] = 0; (1) 
as well as a quantity .E(energy) that is conserved by it: 

dt du a 

For example, if there is a sympletic form and hamiltonian H such that 

V a = {H, u a } 

for the induced Poisson bracket, the invariant measure is the PfafHan (square 
root of the determinant) of u> : 

H = Pf LU. 

Then ([TJ is just Liouville's theorem. Energy will be conserved if 

{H,E} = 

This is automatic if the hamiltonian itself is the energy; but this is not necessary. 
2.2 Dissipation 

A more realistic description could include a dissipative force. We will assume 
that the dissipative force is a gradient of energy; i.e., 



OUh 



dE 
ib 

for some positiv^J 1 ^. The resulting equations 

du„ 



r a + v a 

at 



imply that energy decreases: 



dE _ dE dE < Q 
dt a du a dujj ~ 



2.3 Fluctuation 



We assume a standard model of fluctuating force: a Gaussian r] a (t) with zero 
correlation time ( "white noise") 

<Va(t)Vb(t / ) >= 2D ab 5{t-t'). 

The correlation tensor Z? a6 must be positive as well. The time evolution is 
now a system of stochastic ordinary differential equations 



x That is, r af) is symmetric and T ab u a u h > for all u a . 
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dUg 



Ta + Va+Va- 



The probability density P(u, t) for the random variable u will then satisfy 
a partial differential equation: a diffusion equation with drift. The total proba- 
bility 



P[idu = 1 



so that d ^gP must be a total derivative. This leads us to the Fokker-Planck- 



Kramers equation 



dP _ d 
dt du a 



dP 

(i{D ab —-(r a + V a )P 



2.4 The Boltzmann Solution 

The function 



P = e~P E 



(2) 



is a solution to this equation provided that the dissipation and fluctuation 
tensors are proportional (the Einstein-Smoluchowsky relation): 



For, in this case 



Moreover, 



f3P>ab — r a 6 



D ab ^ = r a p 

ou b 



(3) 



d\uV a P] d\uV a ] dP 

since the ideal evolution V a preserves the measure n and conserves energy. 
If, in addition to (|3j the convergence condition 



Z = j e- f)E [idu < oo (4) 



is satisfied, we can normalize the solution and get equilibrium probability dis- 
tribution 



Z 

If ([3} or Q is violated, the equilibrium is not given by the Boltzmann 
distribution. We will exhibit an example below where ([3} holds but not (|3}. In 
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our case, there is another static, normalizable solution which we will determine 
explicitly. 

In general, there may be no equilibrium state; or the system might reach a 
steady state which dissipates energy at some constant rate. 



2.5 The Adjoint Equation 

Suppose that the Einstein relation holds, so that e~^ E is a static solution. If we 
make the change of variables 



P = e-? E 



Q 



the FPK equation becomes 



dQ 
dt 



1 d 
fi du a 



nD ab 



dQ_ 

du b 



[T a - V a ] 



dQ 
du a 



(5) 



Even when the Boltzmann solution is not normalizable, this is a good starting 
point to search for a normalizable static solution. 



2.6 Fast Dynamics 

If the dynamics V a is much faster than the dissipation and fluctuation effects, 
the details of the vector field do not matter: the system will wander around in 
phase space and fill it. A vestige of the ideal dynamics survives: the invariant 
measure /i. (A kind of micro-canonical ensemble.) In this limit we get the 
Fokker-Planck equation 

BP 

^ = d a [p {D ab d b P - T a P}] . 
Again, if the Einstein relation holds, we have the adjoint equation, 
dQ _ 1 d 

dt jJL du a 



dQ_ 

du b 



. dQ_ 
du a 



3 Stochastic Geodesic Motion on Groups 

There are several examples in physics where we are interested in geodesic motion 
on Lie groups. Arnold's observation that the Euler equations of an ideal fluid 
are geodesic equations on the group of incompressible diffeomorphisms is the 
most important example. Navier-Stokes equations follow from adding a dissi- 
pation and many standard discussions of turbulence involve adding a Gaussian 
random force[4] . Self-erasing random walks can reformulated as random walks 
on the diffeomorphism group of the circle, the Schramm-Loewner Equation j5]- 
This dynamics takes place on a generalization of the upper half-plane, the space 
of univalent functions. A much more elementary example would be rigid body 
motion; random motion of rigid bodies have been used to model dust grains 
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in astronomy [6]. We will now give a general framework for this theory,. The 
ingredients are the structure constants f£ b of the Lie algebra and three sym- 
metric positive tensors G ab , T ab , /^characterizing the energy, dissipation and 
fluctuations respectively. The main result is that the Boltzmann distribution is 
the equilibrium distribution only for unimodular Lie algebras: the trace of the 
matrices in the adjoint representation must vanish. 

We will study in detail the simplest example that is not unimodular. There 
is an equilibrium distribution, but it is not the Boltzmann distribution. 

3.1 Geodesic Motion as Hamiltonian Dynamics 

Geodesies on a Riemann manifold form a hamiltonian system (7] . If the mani- 
fold is a Lie group, and the metric is invariant under the left action of the group 
on itself, the evolution of the tangent vector ("momentum") is independent of 
the position and can be studied separately |8, !) . The Poisson Brackets of veloc- 
ity components are just the commutation relations of the Lie algebra in some 
basis: 

{Va,V b } = fa b V c . 

A left-invariant metric on the group is the same thing as a quadratic form 
on the Lie algebra (":energy") 

E = l -G ab v a v b . 

The geodesic equations follow by the usual rules of hamiltonian mechanics: 

^ = {v a ,E} = f ab G bd v c v d = V a . 

In the special case of the rigid body, the Lie algebra is SO(3); the tensor 
G ab is the inverse of moment of inertia and the geodesic equations are the Euler 
equations of the rigid body. If the metric on the group is also right invariant, 
the tensor G ab is isotropic 

r ab G bd + f d ab G bc = {). 

and the momentum is a constant. This is true of the isotropic rigid body 
(equal moment of inertia in all directions) but is not usually the interesting case. 

More generally, a function H(p) where p 2 = G ab v a Vb can be chosen as the 
hamiltonian. The equations of motion become 

— = {v a ,H} = — x -^f c ab G bd v c v d = V a 
at p 

The solutions are still geodesies, differing only by a constant reparametriza- 
tion: the time variable gets mutliplied by a constant (which can depend on 
energy). We will see that an unconventional choice H(p) = — - will simplify the 
FPK equation in the special case we study in detail below. 
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3.1.1 Invariant Measure 

The obvious measure of integration dv = dv\ A dv^ • • • is only invariant if 

dV a 



For a Hamiltonian H, 

dH 

dv b 



V a = {v a ,H} = f c ab v. 
so that 



dVa =a dH 

dv a Jab dv b - 

A Lie algebra is said to be unimodular if the trace of the structure constants 
is zero: 



r ab = o. 

This is the condition for the measure dv on the Lie algebra to be invariant 
under the adjoint action. The Haar measure on the corresponding Lie group 
will be both left and right invariant. Any semi-simple Lie algebra is unimodu- 
lar; such as the rotations SO(n), linear canonical transformations Sp(n), unitary 
transformations U(n) or any products of these groups. 

An example that is not unimodular is the Lie algebra of the affine group: the 
only non-abelian Lie algebra in two dimensions. Typical examples are solvable 
and nilpotent algebras or Lie algebras containing them as subalgebras. When 
the algebra is not unimodular, there is still a measure fidv that is invariant 
under the adjoint action: it is not anymore just dv. 

3.2 Dissipation 

Adding a fluctuation that is a gradient of energy leads to an Ohmic force 

T a = -T ad G db v b = T b a v b 

and 

— = -r> 6 + f c ab G bd v cVd . 

3.2.1 The Navier-Stokes Equation 

A special case is the Navier-Stokes equation. The divergence free vector fields 
form a Lie algebra under the usual commutator. The energy is just the L 2 -norm 



If 1 

E = - I ViVidx = -G ab v a v b . 
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The tensor G ab is just the Dirac delta function. 

The dissipation tensor, thought of as a quadratic form on vector fields is the 
iJinorm: 



T a v a v b = J diVjdiVjdx 
This leads (after some calculations [9]) to the Navier-Stokes equations 

d 2 v l - VjdjVi - dip 



dvi _ 

— — a Vi 'uj i 
where the pressure pis determined from the constraint 



d t Vi = 0. 

Replacing the space of incompressible vector fields by a finite dimensional 
Lie algebra allows to study simpler models of this important physical system. 

3.3 Fluctuation 

Adding a random force leads to a Langevin equation 



dVg 

dt 



-T° a v b + f c ab G cd v b v d + Va 



The Fokker-Planck-Kramers equation becomes 



dP _ 1 d 
dt ii du a 
There is a static solution 



if 



M { D "bg^+ ( T a v >> ~ fa b G cd v b v d ) P 



Pb(v) 



l3D ab = T ab . 

This can be interpreted as a probability distribution if 



/ 



[ie & 2 dv 



converges. 
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3.3.1 The Randomly Forced Navier-Stokes Equation 

In fluid mechanics, the Navier-Stokes equation with random forcing 
du ■ 

= d 2 Vi - VjdjVi - dtp + rji(x, t) 

is often used to model turbulence. The correlation tensor of fluctuation is 
chosen to be translation invariant 



b{p)d P . 



The function D{p) = p 2 would satisfy the condition (|3j. But it is not the one 
usually used in the literature on turbulence. Generally speaking, dissipation is 
important at short distance scales (viscous force is proportional to the gradient 
of velocity) while fluctuations are believed to pump energy into the system at 
large distance scales. Thus D(p) — |p|~ Q , with a negative power of momentum 
is more reasonable^] . So an equilibrium of the Boltzmann type is unlikely to 
exist in fluid mechanics: instead, a steady state solution in which energy is 
dissipated at a constant rate is more likely to be the correct answer. The theory 
of Kolmogorov that leads to scale invariant velocity correlations, in one such 
example. 



4 The Langevin Equation on a Hyperboloid 

It is often useful to replace an infinite dimensional physical system by a finite 
dimensional toy model, which still has some of the basic structure intact. In our 
case, it is important that the underlying Lie algebra be non-abelian. Otherwise 
it simply reduces to the standard Brownian motion in Euclidean space. Thus 
to get something non-trivial, the Lie algebra must be at least two dimensional. 
The only non-abelian Lie algebra in two dimensions is 

{vo,v±} = Vi. (6) 
The corresponding Lie group is the set of triangular matrices 

f a a\ \ 

low 

with ao > 0. The upper half-plane 

U = {(ao,Oi)|ao > 0} 
parametrizes such matrices. The group multiplication law is 

(a ,ai)(a' ,a' 1 ) = (a a' , a a[ + a x ). 

This group is often called the affine group, as it acts on the real line by 
translations and scaling (affine transformation) 
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(a ,ai) = a t + a\. 

The natural geometry on the upper half plane is the Poincare' metric, which 
has constant negative curvature: the simplest geometry with negative curva- 
ture. There is a one-one correspondence of the hyperboloid with the upper 
half-plane;the induced metric on the hyperboloid is just the the Poincare' met- 
ric on the upper half-plane. A particle constrained to move on the hyperboloid, 
but free of other forces, will move along geodesies. 

The Poincare' metric is invariant under the left action (but not the right 
action) of the affine group. Because it has negative curvature, the geodesies 
do not have constant tangent vectors: only the length is preserved, while the 
direction changes through parallel transport. It is straightfoward to derive the 
geodesic equations [El |9] 

dv 2 

IT = VoVl 

where t is the arc-length. Thee can be thought of as Hamilton's equations 
implied by the Poisson brackets ([6} and the Hamiltonian 

E=\[vl + vl] 
which is just the kinetic energy. 

It is convenient to use a kind of hyperbolic analogue of the polar co-ordinate 
system in velocity space 

v — — ptanh6>, v\ — pesech# (7) 

in terms of which 

<W 

di = P 

d -f = o. 

dt 

The "angular" variable 9 has the range — oo < 9 < oo.lt is canonically 
conjugate to the radial (or "action") variable p 

{P,6} = 1 

under the above Poisson bracket (Ifjl. The discrete variable e = ±1 is needed 
in addition to cover both signs of V\. The solutions are semi-circles in the 
(vq, vi)plane: 

vq = — ptanh pt 
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v\ = pesech pt. 



Such a geodesic motion on a space of negative curvature models the Euler 
equation of an ideal fluid, which are geodesic equations on the group of volume 
preserving diffeomorphisms. Indeed, the afHne Lie algebra is a subalgebra of the 
incompressible vector fields. In a plane for example, the pair of vector fields 



d d d 

—x 

span such a subalgebra 



X dx ^ dy ' ^ dx 



5 Dissipation 



We now modify the above system to include a dissipation, to model the Navier- 
Stokes equation 



dvg 
dt 



-jv - Vx 



dvi 
~dt 



-7W1 + v a v\ 



In terms of hyperbolic polar co-ordinates p, 9 as above ([7f, 



dp 
dt 



-IP, 



d6 
dt 



The energy E(p) — ^p is then monotonically decreasing: 



dE 
~dt 



-2jE 



But the equations can still be solved analytically. The solution that passes 
through the point with hyperbolic polar co-ordinates at time t\ is, 



v = -pie 7(4 *^ tanh 



Vi = piexe 7 ^ *^sech 



Pi- 



h + Pr 



1 _ e -7(*-*l) 

7 

1 _ e -7(*-*i) 

7 



As t — * 00 , the velocities tend to zero; as t — > —00, v\(t) — > and Vo(t) 
—00. For intermediate values it roughly follows the semi-circle. 
If we choose as hamiltonian H(p) = — -,we would get instead 



dp 
~dt 



d9 k 
-IP, -jl = H KP) = -o 
P 2 



dt 



which is just as easily solvable. 
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6 Random Forcing 

Adding a random force with 

< *?a(t) >= 0, < J7a(t)7?b(t0 >= 2D<5 ab 5(t - t'). 
we get the Langevin equation on a hyperboloid 

— = - 7Uo - Vl + m (8) 

dvi 

— = -jvi + v a vi + 772 
at 

As noted earlier, the affine Lie algebra is not unimodular so the invariant 
measure is not dv. It is instead dpd8 = dv v i V = pdv^dvi . The vector field V a 
arising from a hamiltonan H and Poisson Brackets ([6} preserve this measure: 

9H OH d . „. n 

V = Vi = w x — — = 0. 

OVi OVq ov a 

The above (j^J corresponds to H=|?; 2 . The choice H = ~ yields equivalent 
geodesic equations, and is more convenient for solving the FKP equation. 
We are led to the Fokker-Plank equation 

^ = ^d a \n {Dd a P + (jv a - V a )P}} . 

The function e~' ^~ is a static solution if 

0D = r 



But the integral 



is logarithmically divergent near v 1 — 0. In canonical co-ordinates the mea- 
sure is just the Liouville measure 

[id 2 v = dpdO 

In this language, we have a linear divergence in the angular co-ordinate: 



dJS 



Thus the Boltzmann distribution is not the correct equilibrium solution. 
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6.1 The Adjoint Equation 

After the change of variables 

P = e~^Q 
the FPK equation becomes its adjoint 



dQ 

at 



D d 

H dv a 



[i 



dQ 

dv b 



+ [~lV a - V a ] 



5Q 

dv a 



(9) 



Notice that the the first term on the r.h.s. is a kind of Laplacian, but is not 
the Lapalce-Beltrami operator: the symplectic volume element implied by the 
Poisson brackets appears in place of the Riemannian volume element . In polar 
co-ordinates 



dQ 
dt 



D 



d 2 Q 
dp 2 



dQ Dd_ 
lp ~dp~ + p^rn 



We want a solution for which 



dpe 



cosh 9— 



dOQ(p,0,t) 



H'(P) 



dQ 

do 



converges. In particular, it cannot be a constant in the ^variable. 

If the hamiltonian is H = — - the last term would be \ §j? and the equation 



would be solvable by separation of variables: 

Q(p,0,t) = Qi(0)Q2(p)e 



-Dfot 



d_ 
80 



cosh 



} dQi 
d9 



8Qi 
l ^0~ 



PiQi = o 



d 2 Q 2 a dQ 2 13 



^Q 2 + /3 2 Q2 = 



dp 2 dp p* 

where a = ^ . The separation constants /3i , /?2 are eigenvalues of these dif- 
ferential operators; they depend on a as well as on discrete 'quantum numbers' 
labelling eigenfunctions. For a static solution /3 2 = 0. 



6.2 Angular Equation 

Change to 



u = tanh 



d 2 Qi , „ dQ 1 



du 2 



+ 2a- 



du 



+ 



1 — u- 



;Ql = 



13 



Qi = e~ au y/l - u 2 (j)(u) 



(1 - u 2 )(j)" - 2u(j)' + 



01 -a 2 (l - u 2 ) 



This is the Prolate Angular Spheroidal Wave EquationjlO]. 

A simple solution is <j>(u) = yl — u 2 when a = 0,01 = 2. Then Qi[u] = 
(1 — u 2 ) — sech 2 9 is integrable. This belongs to a different branch from the 
non-integrable constant solution, for a — 0. When a ^ 0,this solution continues 
to a solution Sii((9)integrable in 9, with (3i > 2. 



6.3 Radial Wave Function 

For a static solution, 

d 2 Q 2 a dQ 



op z op p z 

The solution is a confluent hypergeometric function. The solution is 



Q2 



l + y/l+4fji 



F 



1 + VI + 4/3; 



If /3i = 0,a solution is constant, which is not normalizable in the angular 
co-ordinate. 

Now suppose a = 0, Pi = 2,for which the angular integral converges. Then 
there is a solution that is normalizable in the radial variable as well. 



1 



1 . 
-e 4 eri 



rpp 



fp P 



Thus, although the Boltzmann distribution is not a normalizable solution, 
there is another solution that is static and normalizable. But it breaks rotation 
invariance spontaneously!. This solution can be continued to 0i > 2, and re- 
mains integrable, when the hamiltonian is not zero. Thus rotation invariance is 
spontaneously broken and moreover, the most likely value of energy is not zero. 

6.4 Equilibrium Solution 

We have putting all of the above together, in the limit of k — when the 
hamiltonian has a small effect, 



P(p, 9)dpd9=e-% p2 



or in the original co-ordinates 



1 , 

— e 4 erf 



2 



sech 9dpd9 
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P(vo,Vi) 



dvodvi vi _ti 



This is plotted in the figure. 



1 . 

-^^e 4 ert 

ftp 



ftp 
2 



dv Q dvi 



Hvo,vi) 




When a ^ 0,we still have an explicit solution 



P(vo,vi) 



dvodvi 

Vi 



which has a similar shape. It is easy to check that this is normalizable. 



7 Conclusion 

There is a large literature on Brownian motion in curved manifolds. Mostly, this 
amounts to a study of the heat kernel of the Laplace-Beltrami operator. This is 
justified in the over-damped limit, where velocity (rather than acceleration) is 
proportional to the force. Our work shows that there are unexpected subtleties 
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in the more general case. The velocities do not tend to the expected Maxwell- 
Boltzmann distribution asymptotically, although for short times they might 
appear to do so. 

For fluid mechanics, the phase space is infinite dimensional. It is possible 
to embed the phase space into an infinite dimensional analogue of the upper 
half plane (the space of complex symmetric matrices with positive imaginary 
part), which has a very similar behavior for the velocities. We plan to return 
to this case, which requires much deeper mathematics than used in this paper, 
in a future publication. It is hoped that the equilibrium velocity distribution 
of fluids under random forces is experimentally accessible and gives useful some 
information about strongly turbulent flows. 

8 Appendix: Formulas For Change Of Variables 

v Q = — ptanhO, v\ — pesech^ 
dv a — — dpta,nh9 — psech 2 9d9 
dvi — edpsech9 — pesechO tanh QdO 



dv^ + dvf = dp 2 + /? 2 sech 



2/) j/)2 



9ab 



p 2 sech 2 6» 



ab 



\ cosh 2 9 
p 2 



dvo A dvi = pesechOdp A d9 



A* = 



psech 9 



pdvo A dvi = dp A d8 
Thus these co-ordinates are canonically conjugate: 



{P,0} = 1 
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